Optimum unambiguous identification of d unknown pure qudit states 
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We address the problem of unambiguously identifying the state of a probe qudit with the state of 
one of d reference qudits. The d reference states are assumed pure and linearly independent but we 
have no knowledge of them. The state of the probe qudit is assumed to coincide equally likely with 
either one of the d unknown reference states. We derive the optimum measurement strategy that 
maximizes the success probability of unambiguous identification and find that the optimum strategy 
is a generalized measurement. We give both the measurement operators and the optimum success 
probability explicitly. Technically, the problem we solve amounts to the optimum unambiguous 
discrimination of d known mixed quantum states. 
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I. INTRODUCTION 

It has been shown that N given pure quantum states 
can be unambiguously discriminated provided they are 
linearly independent [11], at the expense of a certain frac- 
tion of inconclusive results if the states are nonorthogo- 
nal. Unambiguous state identification is a variant of the 
pure-state discrimination problem. Here we assume that 
a quantum system of a d-dimensional Hilbert space is pre- 
pared in a definite pure state out of a set of N linearly 
independent states, but that we do not know the states 
in the set. Instead, we are given N reference systems 
each being prepared in another one of the N unknown 
pure states. Our task is to determine the optimum mea- 
surement for unambiguously identifying the state of the 
probe system with the state of one of the reference sys- 
tems, that is to find the particular measurement which 
maximizes the overall probability of a successful identi- 
fication. Despite the complete lack of knowledge about 
the pure states, their probabilistic identification is possi- 
ble due to symmetry properties inherent in the quantum 
mechanical description. 

Up until now state identification of unknown pure 
states has only been considered for the case N = 2. The 
problem has been introduced by Bergou and Hillery 
for the measurement strategy of optimum unambiguous 
identification and was first solved for two unknown qubit 
states, that is for d = 2 0. The identification of two 
unknown pure states has been shown to be equivalent to 
the discrimination of two known mixed quantum states 
d, 01 • The treatment has been extended to the case that 
the relevant states are each encoded into a certain num- 
ber of identical copies of the respective quantum systems 
d, 0, H, H, 0|, and the measurement strategy of state 
identification with minimum error has been derived, as 
well [3, 01 • Moreover, it has been assumed that only 
one out of two possible qubit states is unknown and the 
other is known In addition, a number of schemes 
for implementing the unambiguous identification of two 
unknown qubit states have been proposed [1, HI, [l0[ . En- 
coding of the unknown states in quantum systems with 



Hilbert space dimension d > 2 has also been considered 
for N = 2, treating both minimum error @, [lH and opti- 
mum unambiguous identification strategies of two states 

Sin. 

In this paper we focus on the measurement strategy 
of optimum unambiguous identification. We extend the 
previous investigations to allow for an arbitrary number 
N of linearly independent unknown pure states. Clearly, 
in order to encode these states we need quantum sys- 
tems with dimensionality d > N. Here we shall consider 
the simplest case, assuming that d = N , and that the 
probe qudit is equally likely prepared in any of the d 
reference states. In Sec. II we establish the connection 
between the unambiguous identification of d unknown 
pure qudit states and the unambiguous discrimination of 
d known mixed quantum states. Sec. Ill provides a gen- 
eral representation of the detection operators describing 
the measurement that unambiguously identifies the un- 
known pure states. Based on this result, in Sec. IV the 
optimum structure of the detection operators and thus 
the optimum measurement, realizing the maximum suc- 
cess probability, is obtained. Sec. V concludes the paper 
with a brief discussion of the results. 



II. UNAMBIGUOUS IDENTIFICATION AS 
MIXED-STATE DISCRIMINATION 

The total quantum system we are considering consists 
of one probe qudit, labeled by the index 0, and d refer- 
ence qudits, labeled by the indices 1, . . . d. Let the states 
\i)n denote orthonormal basis vectors for the n-th qudit. 
Then the identity operator in the d-dimensional Hilbert 
space of the n-th qudit is given by 



/„=^|i)„(i|„ (n = 0, 1, . . . ,fi). 



(1) 



The projector Po^„ onto the two-qudit subspace of dimen- 
sion jointly spanned by the eigenstates of the probe 
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qudit and the n-th reference qudit can be written as 

Po.n = h®In = PoT + PS% {n ^ 0), (2) 

where the projectors Pgfj and -Pp^™ refer to the anti- 
symmetric and symmetric part of the joint subspace, re- 
spectively. In the foUowing we shaU sometimes omit the 
symbol denoting the tensor product. From the explicit 
expressions for the antisymmetric and symmetric basis 
states in the two-qudit subspace we obtain the represen- 
tations 



j = l 1=0 



\i)o\i)n - |i)oN)« {i\o{i\n - 01o(j|n 



V2 



V2 



and 
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j=l i=0 



(3) 



(4) 



\i)o\j)n + \j)o\i)n (i|o01n + {j\o{i\n 



V2 



V2 



which show that the ranks of Pq * and Pq^™ are d{d—l) /2 
and d{d + l)/2, respectively. 

Let us now assume that the state of the probe qudit 
coincides with the state of the n-th reference qudit, so 
that the state of the total {d + l)-qudit system can be 
written as 



|7/;„)|Vl)...|V'n)...|^<i). 



(5) 



On the right-hand-side the order of qudits from left to 
right is probe qudit (zcroth position) followed by the 
d reference qudits. If no information is available about 
the reference states other that they are linearly indepen- 
dent, the discrimination strategy must be independent 
of the actual reference states. In order to emphasize this 
state independence, we introduce the density operator 
I ^'n) (^'n I and take its average over the unknown reference 
states. It will uniformly span the symmetric subspace of 
the probe and n-th reference qudit and uniformly span 
the Hilbert spaces of the remaining d — 1 reference sys- 
tems. Prom symmetry considerations it then follows that, 
for the case when the probe matches the n-th reference 
state, the total {d + l)-qudit system is characterized by 
the average density operator 



Pn 



(d+l) 



PoT ' 



(n-l,...,rf), (6) 



in analogy to the treatment introduced for the identifi- 
cation of two unknown states [1, 01 . Here the pre- factor 
arises due to normalization, taking into account the ranks 
of Pq^™ and of the d — 1 identity operators. Obviously, 
since the state of the probe can coincide with any one 
of the d reference states, there are d resulting average 
density operators. 



In order to identify the state of the probe qudit, we 
have to discriminate among the d different density op- 
erators Pn- A measurement suitable to accomplish this 
task can be formally described with the help of the {d+l) 
positive detection operators Hi , . . . lid and 11? which to- 
gether span the identity. 



d 
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Hv = / = /o (g) /i 



Id 



(7) 



The detection operators have the property that Tr(/9„n„) 
is the probability of successfully identifying the density 
operator as p„, while Tr(p„n„i) (m ^ n) describes the 
probability to get an erroneous result and Tr(p„n?) is the 
probability that the measurement result is inconclusive, 
i. e. that the attempt at discrimination fails to give a 
definite answer [l^, [lj| . For unambiguous discrimination 
we have to require that Tr^pnllm) = for m ^ n, and 
from the positivity of the operators p„ and n,„ it follows 
that this requirement can be only met when 

Pnllni = if m ^ n (8) 

[H, . When the d mixed states occur with equal prior 
probability, given by 1/d, the overall success probability 
of the discrimination measurement takes the form 



1 

Psucc = Tr(pnlln) 



(9) 



For determining the optimum measurement, we want to 
find the detection operators n„ that maximize the suc- 
cess probability Psucc under the constraint that the eigen- 
values of the operator 11? = / — X]n=i -^n have to be 
non-negative. 



III. GENERAL STRUCTURE OF THE 
DETECTION OPERATORS 

A. Three unknown qutrit states 

In the following we proceed along the lines of our pre- 
vious investigations on the unambiguous discrimination 
of two mixed quantum states [Tsl . IigI . [TtI and start by 
deriving the general structure of the detection operators 
n„. It is instructive to consider first only three density 
operators pi, p2 and ps given by Eq. ^ with d — 3 
and referring to the unambiguous identification of three 
unknown pure qutrit states. 

To be specific, let us focus on the detection operator Hi 
that unambiguously discriminates the state pi and there- 
fore fulfills the requirement II1P2 = II1P3 = 0. Clearly, 
this requirement is only met when the support of Hi is 
orthogonal to the supports of p2 and ps, or, in other 
words, it belongs simultaneously to the kernel of p2 and 
to the kernel of ps- From Eqs. ^ and (O we can write 
the projectors onto these two kernels as 



Pk, 



pas 
M),2 



Pk, 



(10) 
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Since Ii is common to both, the projector Pqi onto the 
support of Hi can be written as 



(11) 



P/[^ projects onto the subspace spanned by all 



where 

states \(p) that are linear combinations of the eigenstates 
of the operator Pq2^I3, on the one hand, and also linear 
combinations of the eigenstates of the operator Pq | (g) I2, 
on the other hand. We denote the corresponding eigen- 
states by I a;) for Pg"! (g) I3, and by for Pg^ (g) I2, 



2-^/^(K>0b->2-N>2b>0)|fc>3 -> h>, (12) 
2'^/'(K)ob->3-N)3b)o)|fc)2 - \bl), (13) 

with i < j and i, j, fc = 0, 1, 2, where / — 1, . . . , 9 labels 
the nine triples {i,j,k}. For determining \ip) we put 



(14) 



1=1 



where a; and bi are some complex coefficients. Since a; = 
J2i' bi'{ai\h') and 6; = J2i' {hW) , we obtain 



hi=Y. ^i" Y.^bi\ai,){av\hi,.) = hv, Ba„ , (15) 



i"=i 



where Bu = B22 = P33 = 1/2, P12 = P21 = P23 = 
S32 = -1/4, Pi3 = P31 = 1/4, and P„» = if > 4. 
Here, Eqs. and have been used, with Z = 1, Z = 2 
and Z = 3 standing for {z, j, fc} equal to {1, 2, 0}, {0, 2, 1}, 
and {0, 1,2}, respectively. The only non-trivial solution 
of the system of equations given by Eq. (fTSj) with / — 
1,...,9 reads &i = &3 = — c, h2 — c with c being an 
arbitrary constant, while = if / > 4. Upon inserting 
these values into Eq. and applying Eq. we 

arrive at the normalized state 

l^i(3)>= (-1)^ -^(|0)o|l)2|2)3-|0)„|2)2|l)3 

+ |2)o|0)2|l)3- |2)o|l)2|0)3 

+ |l)o|2)2|0)3-|l)o|0)2|2)3), (16) 

which is the single eigenstate of P^^ . In |</3„(d)) the index 
n = 1 refers to the fact that it is the eigenstate of this 
projector and the argument d = 3 to the fact that we 
are dealing with the qutrit case. It should also be noted 
that \ipn{d)) is the completely antisymmetric state of all 
qudits with the one corresponding to the index omitted. 

In a completely analogous way we can represent the 
projectors on the other two detection operators as Pn„ = 
In ® |<y^'n(3))(<y5„(3)| with n = 2, 3, where 

1^2(3))= i-ir 4^(|0)o|l)i|2)3-|0)o|2)i|l)3 
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|2)o|0)i|l)3- 
|l)o|2)i|0)3- 



|2)o|l)i|0)3 
|l>o|0)i|2)3) 



(17) 



1^^3(3))= (-1)3 -^(|0)o|l)i|2)2-|0)o|2)i|l), 

+ |2)o|0)i|l)2-|2)o|l)i|0)2 
+ |l>o|2)i|0)2-|l)o|0)i|2)2). 

Thus we obtain for n = 1,2,3 



(18) 



Pn„ - /„ ® |^„(3))(vp„(3)| = J2 (19) 



fe=0 



(20) 



where we introduced the normalized states 

l4'^) = l^)n®l¥'n(3)) (fc = 0,l,2). 

Since the states IttI'^'') are orthonormal basis states of the 
support of n„, it follows that n„ can be represented as 

2 

n„= ^ a^'^') |4'=))(7rf )| (n= 1,2,3), (21) 



k,k'=0 



where the coefficients al'^'*^ ^ are some complex constants 
that have to guarantee the positivity of the detection 

operators. Eq. ((20)) yields (tt^ 
and 



= Ok,k' (TTm 1 7r„ ') 



) = [^i 



H ) = -T, 



(22) 



for fc = 0, 1, 2, where Eqs. P^ - P^ have been used. 

A note is in place here about the choice of the overall 
signs of the states |(^„(3)). Although they do not af- 
fect the projectors Pn„ , they lead to sign changes in the 

overlaps of the iTr'h states. For any choice of the overall 
signs either one or all three of the overlaps are negative. 
For the purpose of this paper we made the choice given 
by Eq. ^ together with Eqs. (Ull) - (HH]). With this 
all three overlaps are negative in Eq. (|22p. reflecting the 
intrinsic symmetry of the identification problem. 



B. d unknown qudit states 

Now we return to the general task of identifying d un- 
known qudit states or, equivalently, of discriminating the 
d density operators given by Eq. By the same rea- 
soning that, for d = 3, led to Eq. pT|) . the projector onto 
the support of the detection operator n„ is of the form 



PL 



(23) 



Here P/j 



projects onto the subspace spanned by all 
states that can be simultaneously written as linear com- 
binations of the eigenstates of any one of the operators 
Pom ®i where i ^ m^n and m ^ n. In the following 
we shall show that, generalizing Eq. ((TH]), 



Pn 



®|^„(d))(^„(d)|=^K?))(7rW| (24) 



fe=0 
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with 



ki'^) = |fc)n®l<^n(d)) (fc = 0,l,...d-l), (25) 



where in analogy to Eqs. - pB]) 



l<^n(c?)) 



(-1)" 



(26) 



j = 



In Eq. (j26p the sum is taken over all d\ permutations a 
distributing the excitation numbers aj = 0,1, ... ,d — I 
over the system of d qudits that we obtain by omitting 
the n-th reference qudit from the total system of d + 1 
qudits. The latter are written in fixed order, and sgn((T) 
is the sign of the permutation. 

Since (ttI'^'' IttI'^ ') = 5k, k', Eq. ((24)) is equivalent to the 
fact that the detection operators can be written as 



,d), (27) 



where the coefficients aii''^ •* are complex constants sub- 
ject to the constraints n„ > and II? = / — -^n > 0. 
In order to prove Eqs. and ([?7|) . respectively, we first 



use Eq. ^ to calculate p„ 



for m ^ n, yielding 



(28) 



|A:)„^sgn(a)(g)|cT,),Po^_rko)o|cT. 



= 0. 



Here the equality sign holds because the contributions 
of any two terms P'o,'^K>o|j>m and Po,n™b'>oN>m {i ^ j) 
cancel due to the opposite sign of the respective permuta- 
tions. Hence the requirement for unambiguous discrim- 
ination, Eq. ([5]), is met and Eq. ([M)) is sufficient for 
the operator Pn„ to be the projector onto the support 
of the detection operator H„. It remains to be shown 
that Eq. is also necessary. Because of Eq. this 
holds true if Eq. is the only eigenstate of Pn„ i that 
is if P{j^ is an operator of rank one. From Eq. ([6]) it 
follows that for any m with m ^ n we get the represen- 
tation pm ^ In ® p'mi whcre all operators pj„ together 
span a certain Hilbert space /C„. Because of Eq. the 
eigenstates of also lie in /C„, and the requirement 
Pn„Pm = is equivalent to pj„ = 0. This holds for 
any m ^ n. If P^^ is an operator of rank r, this implies 
that to each of the basis states of the Hilbert space /C„ 
jointly spanned by the operators there belong r states 
that are orthogonal to it and lie also in /C„ . Clearly this 
is a contradiction if r > 2. Thus Eq. ([M)) with r = 1 is 
not only sufhcient, but also necessary for Pn„ to be the 
projector onto the support of H„. 

Before proceeding, we derive two important properties 
of the (d-f l)-qudit states IttI'^''). First, Eqs 



and 

show that for arbitrary n the total number of excitations 



in the {d + l)-qudit system described by IttI*^^) is equal 
to fc -I- X]f=o^ Hence for any two states with different 
k the excitation numbers of at least one of the qudits 
composing the total system are different and therefore 
these states are orthogonal, which leads to 

(^L'^kf^>=4,.'(7rW|4'^>- (29) 
Second, while {iT^n''\Tr^n^) — 1, for m ^ n we find that 



(30) 



since the expression 



dl 



E Sgn(o-) sgn (o-' ) (fc I cr„) „ (cr^. I cr^ ) j « | fc) r 



j=o 



is equal to —1/d. Here we took into account that the 
inner products vanish unless a'j — Uj for j ^ m,n and 
CT^j = fc = Gm, reducing the double sum to a single one 
over (d— 1)! permutations which each contribute the value 
1. Eq. (|30p reflects the intrinsic symmetry of our identi- 
fication problem. 



IV. THE OPTIMUM MEASUREMENT 

Now we apply the general representation of the detec- 
tion operators H„, Eq. ([77|) . in order to find the special 

coefficients a'n'^ that determine the optimum operators 
H„ which maximize the overall success probability Psucc, 
given by Eq. ©. Using Eqs. (O and (|77)l a straightfor- 
ward calculation shows that 



Tr(p„H„) = "n'"'^ M'Vnkf ^ 



(31) 



k,k'=0 



d-1 



with 



d-i 



41U^(z|o(fckPoTlfc') 



d+l 



* = 



Sk,k', (32) 



i=0 



where Eq. (j4|) has been apphed in Eq. (f32|) . Thus, we 
arrive at 



d d—1 / d 

n=l k=0 \n=l 



, (33) 



where = al'^''^''. In order to maximize Psucc we need 

(k) 

to determine the largest values a„ which are still in 
accordance with the constraint H? = / — ^ 0, 

•-"^ X^n Hn < I, respectively, where the operators H„ are 
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given by Eq. (|27|) . Since Pgucc does not depend on the On the other hand, from Eq. (|38p we obtain 

non-diagonal elements (Mi'^ "* (fc 7^ ^r maximizing 
Psucc on the given constraint we have to put 



a 



{k,k' 



a. 



k,k' 



The constraint then reduces to 

n=l A:=0 \n=l 



(34) 



(35) 



where from Eq. ((29)) it follows that for different values of 
k the operators within the bracket in Eq. ([35]) act in or- 
thogonal subspaces. Maximizing the success probability 
given by Eq. (|33p therefore amounts to solving d inde- 
pendent maximization problems in the d orthogonal sub- 
spaces belonging to different values of k. Moreover, since 
according to Eq. ([50]) the mutual overlaps (TTm^ Itt^'^^) 
do not depend on k, all these optimization problems are 
mathematically identical and the index k can be omitted. 
Therefore our task is to maximize X]n=i Q;„ = S' with 



arik„)(7r„| = lis < /, 



(36) 



(7r„|7r„) = 1, (7r„j|7r„) = -- (m^n). (37) 

For this purpose we make use of the method developed by 
Chefles and Barnett for the optimum unambiguous 
discrimination of linearly independent symmetric pure 
states. We first observe that d state vectors, with mutual 
overlaps given by Eq. ([57)) . can be represented as 



Vd 



1 '^"^ n 
- ^exp {2111^1^ \ui) 



k„) = ^ |wo) 

1=1 

(38) 

where the states with — Sur denote some 

orthonormal basis in the d-dimensional subspace spanned 
by the set of states {|7r„)}. Indeed, from Eq. (|38p we get 



1 
d 



d+1 



d-l 

E 

(=0 



exp 27ri 



1 1 



(39) 



and Eq. (j37p is immediately recovered using the sum 
rule for a geometric series. With the help of Eq. ([38| it 
is easy to check that the states {|7r„)} belong to the class 
of symmetric states that are covariant with respect to 
the unitary operator U = ^f^Q {'^T^iji) k/)('i*;| £ind 
transform according to U\7rd) — \'^\) and f7|7r„) — |7r„-|_i) 
for n = l,...c?— 1. It has been shown [isl that then 
there exists an optimum operator 11^'' maximizing S and 
possessing the same symmetry with respect to t/, leading 
to ai = . . . = a„ = a. The constraint given by Eq. (I36|) 
thus takes the form 



d 



(40) 



(41) 



n=l 



1=1 



where the relation 'Ylm=\ exp[27ri(Z — Z')n/(i] = d Sw has 
been taken into account. Inserting Eq. ([1T|) into Eq. ([^0]) 
we find immediately that the largest value of a, allowed 
by the constraint, is given by a°^* = d/{d + 1). With 

a(f ' = a°P\ Eqs. 1^ andl^ yield both the max- 

imum success probability for unambiguously identifying 
d unknown pure qudit states, 



popt ]_ 

^'^'^ {d+l)dd-^' 



(42) 



and the optimum detection operators, with IttI'^'') given 
by Eq. 



d-l 



in = l,---,d). (43) 



k=0 



1=0 



Equations (|^^ and represent the main results of 
this paper, generalizing previous results from two qubits 
to d qudits. For d = 2 Eq. ([1^ reproduces the success 
probability 1/6 obtained previously 0, 0| for the opti- 
mum unambiguous identification of two unknown pure 
qubit states. The optimum detection operators 11°^* are 
not projectors and, therefore, the optimum measurement 
strategy is a generalized measurement. 



V. CONCLUSIONS 

In this paper we derived the optimum measurement for 
the unambiguous identification of d unknown pure qudit 
states. Our treatment is based on Eqs. ([^5]) and ([M]) . re- 
sulting from the special structure of the density operators 
to be discriminated, and reduces the optimization prob- 
lem to d independent maximization problems in orthogo- 
nal subspaces of dimension d. The procedure is analogous 
to determining the optimum measurement for the dis- 
crimination of two mixed quantum states in cases where 
the optimization reduces to maximization problems in 
orthogonal two-dimensional subspaces [l^, Il6l Iv^ . It 
should be noted in this context that the correct detection 
operators have previously been derived by Zhang el al. 
using a slightly different method and without explicitly 
working out the optimal success and failure probabilities 

Quantum systems with d > 2 have been proposed as 
carriers of quantum information in various contexts like 
e. g. quantum cryptography [l^, and methods for re- 
alizing general linear transformations on single-photon 
qudits have been theoretically described [2(]||. Moreover, 
the manipulation of biphotonic qutrits [2l| and ququarts 
j22i | has been experimentally demonstrated. Apart from 
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possible applications in quantum information, our inves- 
tigations are also of interest with respect to the theory of 
optimum unambiguous discrimination of more than two 
mixed quantum states [1^, 113] , where no examples of ex- 
plicit solutions have been given before. 
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